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We present analysis of a single channel interacting quantum wire problem in the presence of 
spin-orbit interaction. The spin-orbit coupling breaks the spin-rotational symmetry from SU(2) 
to U(l) and breaks inversion symmetry. The low-energy theory is then a two band model with a 
difference of Fermi velocities 5v. Using bosonization and a two-loop renormalization group procedure 
we show that electron-electron interactions can open a gap in the spin sector of the theory when 
the interaction strength U is smaller than Sv in appropriate units. For repulsive interactions, the 
resulting strong coupling phase is of the spin-density-wave type. We show that this phase has 
peculiar emergent topological properties. The gapped spin sector behaves as a topological insulator, 
with zero-energy edge modes with fractional spin. On the other hand, the charge sector remains 
critical, meaning the entire system is metallic. However, this bulk electron liquid as a whole exhibits 
properties commonly associated with the one-dimensional edge states of two-dimensional spin-Hall 
insulators, in particular, the conduction of 2e^/h is robust against nonmagnetic impurities. 


I. INTRODUCTION 

The discovery of the quantum spin Hall insulatorflUS 
sparked the realization that gapped phases of matter with 
identical bulk spectra are not all equivalent. Indeed, if 
an insulator can be characterized by a nonzero topolog¬ 
ical invariant, it hosts gapless modes at its edge which 
are robust to perturbations that respect the anti-unitary 
symmetries of the system. This topological protection 
crucially depends on a finite gap in the bulk. While 
the topological classification of non-interacting systems 
is well established,!^ the role of strong interactions is still 
a matter of ongoing research. 

A prime example of a strongly correlated system are 
electrons in a quantum wire which ar e a g ood realization 
of a Tomonaga-Luttinger liquid (LL).!^The distinctive 
feature of this state of matter is that the elementary ex¬ 
citations have no relation to free electrons, but rather 
are described by collective plasmon modes. Addition¬ 
ally, these plasmon modes carry spin and charge inde¬ 
pendently, a phenomenon known as spin-charge separa¬ 
tion. As a consequence of the collective nature of the ele¬ 
mentary excitations, even weak interactions between the 
electrons have profound consequences for the quantum 
state of the system. If interactions become strong, they 
can lead to a strong coupling regime where spectral gaps 
are generated dynamically without spontaneous breaking 
of any continuous symmetry. 

A particularly fascinating example of such a dynam¬ 
ically generated state is the Luther-Emery liquid^, in 
which the charge sector remains critical (gapless), how¬ 
ever, the spin degrees of freedom acquire a gap. This 
can quite naturally occur for attractive interactions, in 
which case the pairs of spins form singlets and t he sy s- 
tem exhibits many properties of superconductivity.*^*^ If 
the interactions are repulsive, however, the spin-gap may 


only form if SU(2) spi n-ro tational symmetry is somehow 
broken in the system.*^ While this case has the same 
thermodynamic spectrum as the spin gap that appears 
for attractive interactions, the states are rather differ¬ 
ent, with the dominant correlations in the repulsive case 
being of the spin-density- wave (SDW) type. 

In this work, we will argue that the spin-gapped sys¬ 
tem with repulsive interactions is topologically nontriv¬ 
ial. This topological state is peculiar in the sense that 
(i) it emerges as a strong coupling phase of the origi¬ 
nal model and (ii) even while one sector of the theory 
is gapped the other remains gapless so that the overall 
electron liquid is not in a gapped phase. Nonetheless,the 
system exhibits properties usually associated with topo¬ 
logical insulators. First, zero-energy edge modes with 
fractional spin emerge at the boundary of a finite system. 
We note that this property has recently been predicted in 
a related one-dimensional model with spin-triplet pairing 
in Ref. [101 while various topological properties of gapless 
states have been discussed more generally in Refs.iTTIand 
im Second, since electrons carry both charge and spin 
the gap in the spin sector affects the whole electron liquid 
and leads to unusual transport properties. In particular, 
we find that the bulk of the wire is robust against non¬ 
magnetic impurities as long as interactions are not too 
strong [for Kc > 3/4]. We note that these transport 
properties are inherent to the SDW phase and have also 
been discussed in Ref. M and [m Throughout the pa¬ 
per we will compare the transport properties of the two 
realizations in more detail. 

The study of strongly interacting one-dimensional sys¬ 
tems is by no means academic. Over the last several 
years there has been remarkable progress in nanotechnol¬ 
ogy which has led to an explosive growth of experimen¬ 
tal work on low-dimensional systems. Single and mul¬ 
tichannel one-dimensional conductors or quantum wires 
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can now be manufactured in a controlled fashionfi^H^ 
These systems provide a fertile ground for laboratory ex¬ 
periments of strongly interacting systems. 

However, as interactions in quantum wires are natu¬ 
rally repulsive, the SDW phase can only be realized if 
the spin SU(2) symmetry in the systems is broken. In 
a realistic setup, this is naturally achieved by the pres¬ 
ence of spin-orbital interaction. In experiment, ballis¬ 
tic quantum wires are created in a two-dimensional elec¬ 
tron gas (2DEG) by cleaved edge overgrowth. They are 
therefore naturally subject to spin-orbit coupling (SOC) 
which breaks the SU(2) symmetry. The SOC arises due 
to the asymmetry associated with the potential that con¬ 
stricts electrons to the two-dimensional plane, the so- 
called Rashba SOC.^^The asymmetry and therefore the 
Rashba SOC can be further controlled by applying ex¬ 
ternal gate voltage.^^Mlll addition to the noted asym¬ 
metry due to the confining potential (which include the 
quantum-well potential that confines the electrons to the 
2D layer as well as the in-plane potential that forms 
the quantum wir^^ spin-orbit interaction is inherent 
to semiconductors of zinc-blende or wurtzite structure 
lacking a center of inversion. This leads to the so-called 
Dresselhaus SOC.^In this paper, we consider the situa¬ 
tion where the Rashba term is tuned to be much stronger 
than the Dresselhaus term. This limit can be achieved 
experimentally by applying a sufficiently strong backgate 
voltage. 

The effect of SOC on interacting one-dimensional 
systems has bee n discussed extensively in the 
literature.^i^^i^^2lII3l2il However, so far it was be¬ 
lieved that the SDW state can only be realized with 
the help of an additionally applied magnetic field which 
explici tly b reaks the time-reversal symmetry of the 
systen^^^ or by fine tuning a modulated Rashba 
SOC.™l 

In contrast, we show in this work that the SDW state 
can be realized in quantum wires with strong spin-orbital 
interaction even without additional perturbations such as 
magnetic fields as long as the spin-orbit coupling leads to 
different velocitiepSElMl Fermi points of the low- 

energy bands. In this case, we find the criterion for the 
gap opening that the dimensionless velocity difference Sv 
should be larger than the dimensionless Hubbard U in¬ 
teraction. This constitutes one of the main results of our 
work; throughout the paper we will critically contrast 
our results with previous work on spin-orbit coupling in 
quantum wires in order to explain how we obtained a 
different answer. Having shown that the SDW may be 
realized, we will then discuss the topological properties 
of this state. In particular, we will show that zero-energy 
edge modes emerge at the boundary of a finite wire, and 
electrical conduction in this state is insensitive to impuri¬ 
ties as long as time-reversal symmetry remains unbroken. 
This constitutes the second main result of this paper. 

The paper is structured as follows. In Sec. |n) we first 
use a specific tight binding model to study the effect of 
spin-orbit interaction on the spectrum of non-interacting 


fermions. After having established some basic under¬ 
standing, we summarize results of previous studies and 
contrast them to this work. In Sec. m we use the re¬ 
sults from the noninteracting case to formulate an ef¬ 
fective low-energy theory in the presence of interactions. 
We proceed to bosonize this model which takes interac¬ 
tions into account exactly and discuss the influence of 
SOC by deriving the renormalization group (RC) flow in 
Sec. |IIIB[ We find that the RC flows to strong coupling 
in a certain parameter regime which leads to the opening 
of a spin gap. The nature of the strong coupling fixed 
point is discussed in Sec. [Ill C| and we study the topolog¬ 
ical properties as well as the effect of disorder in Sec. |W 
A summary of the present work is presented in Sec. V 
which is followed by two appendixes. In Appendix X 
we present the bosonization and refermionization con¬ 
ventions used throughout the work, while Appendix 
presents details of the renormalization group procedure. 
Throughout the paper, we use units where h= 1. 


II. TIGHT BINDING MODEL FOR ONE 
DIMENSIONAL ELEGTRONS IN THE 
PRESENGE OF SPIN-ORBIT-GOUPLING 

We consider spinful fermions confined to one spatial 
dimension and subject to a Rashba spin-orbit-interaction 
at incommensurate filling. The Hamiltonian of the model 
is 

H = Ho + Hso + Hint ■ (1) 

The bare hopping Hamiltonian is given by 
Ho=-tY^ + H.c. 

+ H-C. - ^iN 


( 2 ) 




= - 51 cos(2fc) -b ^J\ cj, . 

k,a 


Here Cj^a- destroys an electron with spin a at site 

j. The hopping amplitude is denoted by t for nearest 
neighbor and t' for next-nearest-neighbor hopping. We 
use dimensions where the lattice spacing oq = 1 and we 
assume periodic boundary conditions. 

The Rashba-type SOC reads as 


Hso = - 51 + H-C- 


= 2a^sin(A:)cJ,_^ 


(3) 




The Rashba SOC with coupling strength a breaks the 
SU(2) spin-rotational symmetry down to U(I) but pre¬ 
serves time-reversal symmetry. Here ct*, with i G 
{x, y, z}, denotes the set of Pauli matrices in spin space. 
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Finally, the Hubbard interaction is given by 

Hint = UY^ • ( 4 ) 

3 

with the coupling constant U and the electron density 
operator 

Throughout this work, we assume the hopping am¬ 
plitudes as positive, t,t' > 0 and repulsive interactions, 
U >0. 


A. Effects of spin orbit coupling on the spectrum 
of noninteracting electrons 


First, we discuss the effect of SOC on the spectrum 
of noninteracting fermions without next-nearest-neighbor 
coupling i.e we set t' = 0. The non-interacting part of 
the Hamiltonian Q then reads as 

= [-2^ cos(fc) -k 2acr sm{k)] Cfc ,,, (5) 

k,(j 

where cr = ±1 are the eigenvalues of Using the har¬ 
monic addition theorem, this can be recast into the form 

Ho=- 2i'^cos{k - <yqo)cl ,^Ck,a ■ ( 6 ) 

k,(T 


Here, i = is the renormalized hopping ampli¬ 

tude and go = arctan(a/t). 

We conclude that in the absence of next-nearest- 
neighbor hopping SOC renormalizes the hopping ampli¬ 
tude and shifts the spectrum by a constant momentum 
±goj for spin up or down respectively. However, this shift 
can always be removed by a spin dependent gauge trans¬ 
formation and therefore has no observable effect on the 
thermodynamic properties of the system [they depend 
only on the spectrum]. 

This statement can be made explicit by considering the 
Hamiltonian in Eq. (§ back in real space: 

i/o = - 2<E = -2i E ’ (7) 

j.o- 3,cr 


where we defined the fermion operator 
The model in the presence of SOC is therefore unitarily 
equivalent to a model without SOC but with renormal¬ 
ized hopping parameter.!^ Note that this statement re¬ 
mains true in the presence of interactions since the trans¬ 
formation leaves the density rij^cr = c] o-O.o- = 
invariant and thus does not change the form of the inter¬ 
action term 'Hint- For repulsive electron-electron inter¬ 
action, the system therefore would be in the Luttinger 
liquid phase (see Sec. HIB), as it would be in the ab¬ 
sence of SOC. 

This possibility of gauging out the Rashba term was 
presented recently by Goth and Assaad in Ref. |34| How¬ 
ever, it is curious to note that this Hamiltonian is iden¬ 
tical to t hat of two spinless fermionic chains in a mag¬ 
netic fieldThe Rashba term in the former model is 


equivalent to the orbital magnetic field in the latter; the 
inter-chain coupling in the ladder model being equivalent 
to a Zeeman term in the Rashba case, were such a term 
to be present. 

In the ladder realization of the Hamiltonian, it is fairly 
clear why the magnetic field can be gauged out unless 
the chains are coupled: a pure one-dimensional object 
can have no orbital motion, and therefore can not feel 
the orbital effects of a magnetic field. In its essence, the 
reason why spin-orbit may be gauged out in the Hubbard 
chain is analogous: There is no ’’orbital” motion possible, 
so the spin-orbit may only couple as a pure gauge [one 
should be careful, however, that this is not to say there 
are no observable effects of spin-orbit (see Ref. jSllfor the 
analysis)]. 

However, as any physical manifestation of a one¬ 
dimensional wire is necessarily due to a confining poten¬ 
tial in the transverse directions, the Rashba spin-orbit in¬ 
teraction may nevertheless couple to the wave functions 
in these transverse directions to affect the system in a 
thermodynamic manner. This was a nalyzed 15 years ago 
in a series of papers by Moroz et alplElIIll crux of 
this work is that non-trivial effects occur when the parity 
(inversion) symmetry of the wire is broken in an essen¬ 
tial way (assuming time-reversal symmetry is preserved). 
This broken inversion symmetry may be intrinsic to the 
wire itself, or due to the confining potentials. While we 
refer to these original papers for an overview of the mate¬ 
rials physics aspect of this, it turns out one may capture 
this behavior with a simple toy model, Q, with next- 
nearest neighbor hopping. 

When t' 7 ^ 0, the spectrum of Q is given by 

f-crik) = —2icos{k — aqo) — 2t'cos(2fc). (8) 


We find two bands characterized by the z-component 
of the spin that are shifted by an constant momentum 
go- The band structure is depicted in Fig. Notice 
that the inversion symmetry is broken due to the SOC 
ecr{k) 7 ^ €a{—k). While this is also true when = 0, in 
the presence of next-nearest neighbor hopping this leads 
to different Fermi velocities at the Fermi points of each 
band [cf. Fig. 2], and thus the symmetry can not be 
restored by a trivial gauge transformation. Meanwhile, 
time-reversal symmetry of the model is still preserved 
and therefore eo-(fc) = e-a-{—k). 

We can find an analytical estimate for the velocity dif¬ 
ference in the limit when and the chemical poten¬ 

tial is tuned to the bottom of the band. In this case, we 
can expand the spectrum around k = aqo. We hnd the 
Fermi points of the band a determined by the equation 
£a-{ki^ 2 ) = M and the Fermi velocities: 


_ dcajk) 
dk , . 

k—ki^2 

= ± 2^J(t'Y sin^(2go) + + cos(2go)) 

-I- 4t'crsin(2go). 


(9) 
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FIG. 1. Spectrum in Eq. ^ of the hopping model in Eq. Q 
in the noninteracting case. The bands are labeled by the 2 - 
component of the spin and shifted by the constant momentum 
go- The energy in units of the hopping parameter t is shown 
for the parameters a/t = 0.3 and t'/t = 0.5. 


where fl = 2t' cos{2qo)+2i+fj.. Due to the preserved time- 
reversal symmetry, the Fermi points and Fermi velocities 
are not independent but rather = fci, etc. 

The velocity difference at the Fermi level is given by 




(t'AT 


( 10 ) 


Notice that the velocity difference vanishes either in the 
absence of SOC, a = 0, or without next-nearest-neighbor 
hopping, t' = 0. For weak nearest-neighbor hopping and 
weak SOC, it is of the order of Sv ~ at'/t^. 

To summarize, by using an explicit hopping model for 
one-dimensional fermions in the presence of SOC, we 
have identified two main effects of SOC on the spectrum 
of noninteracting electrons. First, it breaks the SU(2) 
spin-rotational symmetry and therefore lifts the spin de¬ 
generacy of the spectrum and second it breaks inversion 
symmetry leading to different Fermi velocities vi ^ V 2 - 
Before we discuss the effect of interaction on the phase 
diagram of the system we now review existing work on 
the topic to put this work in the correct context. 


B. Summary of previous work 

This work considers two interrelated questions: first, 
whether spin-orbit coupling may drive a [otherwise SU(2) 
invariant] single-channel quantum wire with repulsive in¬ 
teractions to a spin-gapped SDW phase; and secondly, 
what are the physical properties of this phase, in partic¬ 
ular, those associated with non-trivial topology. As, par¬ 
ticularly, the former of these questions has a rich history, 
we find it useful to briefly summarize previous results 
relevant to this work. 


The model w e use wa s first described fifteen years ago 
by Moroz et who derived explicitly that in 

general the SOC gives rise to two bands with different 
Fermi velocities. They modeled a quantum wire as a two- 
dimensional electron gas confined in one spatial direction 
by an external potential and derived the band structure 
of the effective one-dimensional model. The spectrum in 
the presence of SOC turns out to have the same form as 
that of our hopping model and they propose the effec¬ 
tive low-energy theory we discuss in Sec. |III| We note, 
however, that it was not until very recentljl^ that it was 
realized how essential this difference in Fermi velocities 
is, as if this effect is neglected, the SOC may be removed 
by a gauge transformation. 

The effect of interactions was also discussed in the 
early works.^ By a simple scaling dimension analysis, 
they concluded that the backscattering interaction term 
potentially responsible for the opening of a spin-gap is 
always irrelevant (for repulsive interactions), and thus 
a spin-gap never opens. However, it turns out that in 
the vicinity of the SU(2) symmetric point, the scaling 
dimension alone is not a good indication of relevance 
as the backscattering interaction at this point is exactly 
marginal.^ One should therefore study this question more 
carefully. 

Five years later, Gritsev et al^ revisited the prob¬ 
lem of the RG of interacting fermions in the presence of 
Rashba SOC. By treating the problem within two-loop 
RG, they concluded that the SDW phase is possible, and 
constructed a phase diagram. However, they start from 
general parameters in the model [meaning interactions 
may explicitly break SU(2) symmetry even in the ab¬ 
sence of spin-orbit coupling], and it is very difficult from 
their work to determine the line in parameter space the 
physical situation where the SU(2) is broken only by spin- 
orbit. In fact, as they do not consider different velocities 
in the two bands, it must be true within their model that 
the spin-gap state is never realized on this line. 

A few years later, Schulz et al. revisited this problem 
within one-loop RG, but treating the marginal backscat¬ 
tering interaction much more carefully.^ Through this 
analysis, they obtained a Berezinskii-Kosterlitz-Thouless 
(BKT) phase diagram (as is the case without spin orbifl^, 
but with renormalized effective parameters, which de¬ 
pend crucially on the velocity difference. Importantly, if 
the velocity difference goes to zero, the bare backscatter¬ 
ing parameters are not renormalized; which is an equiva¬ 
lent way of saying that the SOC may be gauged out. The 
result of this calculation, however, was that the parame¬ 
ters are always renormalized towards the weak-coupling 
side of the phase diagram, implying once more that the 
SDW state is never realized. 

We follow a very similar approach to Ref. [331 but by 
integrating out the gapless charge sector, we treat ex¬ 
actly the forward scattering couplings between the spin 
and charge sectors before applying perturbative RG. Like 
Schulz et al., we also obtain a BKT phase diagram with 
renormalized parameters, but these renormalized param- 
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eters can now be on either side of the phase boundary, 
meaning that under certain conditions (which are derived 
in this work), the spin-gap may appear and the SDW 
phase is realized. 


We note also that there is a lot of work on the realisa¬ 
tion of the SDW state when a magnetic field is also added 
to the system. The authors in Refs. [T3] and [HI consid¬ 
ered the situation when both spin-orbit coupling and a 
magnetic field are present which breaks SU(2) symme¬ 
try completely. They find that the system undergoes a 
phase transition to the SDW state if the magnetic-field 
and spin-orbit axes are orthogonal. Furthermore, if the 
developed spin gap is large enough, the ordering in the 
spin sector can crucially suppress the backscattering of 
electrons of nonmagnetic impurities. On the other hand, 
magnetic impurities can localize the SDW phase and de¬ 
stroy the (near-)perfect conduction properties. This is 
a general property of the SDW phase and will discussed 
in detail in our work in Sec |IV A One key feature of 
the present case is that this phase is robustly protected 
against these magnetic impurities by the time-reversal 
symmetry in the system, or in other words, as long as 
time-reversal symmetry remains unbroken, this state will 
remain metallic. 


Finally, we comment briefly on the topological proper¬ 
ties of the phase we find, which are similar to those of a 
topological insulator (due to the spin-gap), but occur in 
a metallic system, due to the gapless charge modes. Re¬ 
cently, another work appeared by Keselman and Ber^iSl 
looking at exactly this question. Although Keselman and 
Berg looked at attractive interactions, they were inter¬ 
ested in a spin-triplet pairing rather than the conven¬ 
tional spin-singlet pairing, and this spin-triplet pairing 
has exactly the same spin structure as our SDW. The 
topological properties of the two models are therefore al¬ 
most identical. The fractional spin-edge states we find 
are analogous to those in (TUI while we believe our con¬ 
duction properties will also carry over to their model. 


III. INTERACTION EFFECTS 


This yields the density 

n,-, ^ RlR, + LlL, + 4, 

( 12 ) 


where we defined kp = (fci + k’2)l2. While the SOC 
breaks the chiral symmetry of the low- energy model, 
time-reversal is still preserved. For spinful fermions, 
the time-reversal symmetry operation T in real space 
can be represented as T = idyK. Here, Oy acts in 
spin space and K denotes complex conjugation. This 
implies —)■ ^nd Us¬ 

ing the low-energy decomposition in Eq. ( |11[ ) we find 
R^{x) —?► A*(x), L^{x) —?► R*{x) Ri{x) —>■ —L^{x), and 
Li{x) —)■ —R^{x). Using the bosonization conventions 
outlined in Appendix]^ this implies the following trans¬ 
formation properties in the spin-charge basis: 

(Pc(x) -)> (ficix), ^s{x) -^s(.x), 

9c{x) -)> -Oc{x), 9s{x) 9s{x). 


Additionally the Klein factors transforms as —>■ and 

Kl —>■ —K^. 

The low-energy form of the Hamiltonian is given by 


Ro =- ivi J dx (^R\dj:R^ - 

- iv2 J da; (rID^R^ - L\dj;L^'^ , 

Hint =ujdx {R\R^ + l\l^) (rIr^ + 
+U j da; {R\l^lIr^ -b H.c.) . 


(14) 

(15) 


Here, Vi and V 2 are considered as phenomenological pa¬ 
rameters of the low-energy theory which describe the dif¬ 
ferent Fermi velocities at the left and right Fermi point. 
On the basis of our analysis of the hopping model in 
Sec. [H] we expect the velocity difference to be small but 
in principle tunable through the Rashba parameter a. 
The interaction term Hint follows from applying the de¬ 
composition (12) in the Hubbard interaction in Eq. (|^. 


We now construct an effective low-energy Hamiltonian 
based on the previous results on SOC effects. In the pre¬ 
vious section we have learned that the SOC lifts the spin 
degeneracy of the energy bands and leads to two bands 
(tj =t,4.) with asymmetric Fermi velocities Vi ^ V 2 - To 
find the effective low-energy form of the Hamiltonian we 
linearize the spectrum near the Fermi points fci 2 , see 
Fig. H Next, we expand fermionic operators in modes 
that vary slowly on the scale of the inverse Fermi mo¬ 
mentum: 


Cj-^ ^ ^^(x) = R^{x) + ^{x) , 

Od ^ V’t(a;)=H;(a;)e*'==" + L;(x)e-*'=iU 


A. Bosonization 


We study the low-energy theory with the help of the 
bosonization technique using the conventions outlined in 
Appendix The resulting Hamiltonian density of the 
kinetic part reads as 




(7 

6v r 




(16) 


+ y 


where we introduced the difference 5v = Vi — V 2 and 
average vp = ('fi + U2)/2 of Fermi velocities . 
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FIG. 2. Band structure e{k) of the low-energy theory. Due to 
time-reversal invariance the dispersions of the two bands are 
connected as e-f{k) = e_i{—k). The low-energy excitations for 
up spins are right moving particles with velocity wi at Fermi 
momentum fei and left moving particles with velocity —V 2 at 
Fermi momentum —k 2 (analogously for down spins). 


the SOC manifests itself only in the term Hso that breaks 
the spin-charge separation. In the absence of SOC the 
RG flow of the model is described by the well-known BKT 
equations (351 and the corresponding flow is constrained 
by SU(2) symmetry to the separatrix(see Fig.|^. Now, as 
long as the velocity asymmetry Sv is sufficiently small we 
can stay in the spin-charge basis and treat the term "H as 
a small perturbation. Usually, the (marginal) term Tiso 
would be neglected as it only produces small corrections 
under the RG. In the present case we have to keep it since 
the conventional flow is exactly along the separatrix and 
even a marginal term may drive the system into a new 
phase. 

To study the effect of this perturbation on the RG flow 
of the model we will integrate out the (quadratic) charge 
sector and derive the effective action of spin fields. 

First we perform the Legendre transformation from the 
Hamiltonian to the Lagrangian 


C [(fii, d^(pi, dtipi] + H Hi 

= Jdx . 


( 22 ) 


We now introduce the conventional spin and charge 
operators 


= 


n + ‘Pi 

V2 ’ 


= 


n - 


(17) 


In the new basis the Hamiltonian density including 
electron-electron interaction reads as 




2 L 


Kj\l + K-^ 


K.ni + (a,^e) 


9 s 


2(7ra)- 


cos 




■^SO = 


Sv 


2 L 


K — 1-1 _ 

^ + 2TTVF ’ 


(18) 


-t- 

Here, the Luttinger parameters are defined as 

= c,s . (19) 

Due to Galilean invariance of the model it must hold that 


Vf^K^ = VF ■ (20) 

In terms of the original parameters, we have 

9c = -gs = -aoU, (21) 

where we reinstated the lattice constant qq. In particular 
for repulsive interaction U > 0, we find < 1. 

We notice that the charge sector is described by a LL 
state with coupling constants Vc and Kc- The spin sector 
is also a LL but includes a backscattering term that can 
generate a spin gap if it becomes relevant in the RG sense. 
We point out that the different Fermi velocities and thus 


where dt(pi = SH/SUi. Solving the equations for the 
conjugate momentum H yields the expressions 


TT _ 1 c) 

Vp ZVp 

TT _ 1 c, 

Vp ZVp 


(23) 

(24) 


Substituting Eqs. (181 and (24) into Eq. (22) yields the 
resulting Lagrangian density 

, _ 1 r... ,2 9 


2vcKc 

1 

' 2v*K* 




c c 

1 


2vsKs 

9s 


(dtific) - '('c ( 1 - ^ 


(dtipcf - {d^(pcy 


{dt^PsY {d^Ps) 


(25) 


2(7ra)- 

Sv 


cos 




^SC - l^xPc^tPs T ^xPs^tPe] • 

ZVf 

Here we introduced the renormalized coupling constants 
Sv'^' 


K*=KAl- 


dv'i 


V* = vf/K 


(26) 


Next we determine the action in imaginary time 
iS =i j dtdx C[pAxPidtp\ 


t—lT 


—S = Jdrdx ClpAxPydrp] ■ 


( 27 ) 
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The partition function is then 

Sr. = [{dr^cf + {vlf {d^ificY] , 


5., = 


2vp 


9s 


dxdr cos 




2(7ra)^ 
zSv f 

Ssc = ^ J dxdr [d^^PrdriPs + d^^sdr^^c] ■ 


(28) 

(29) 

(30) 

(31) 


Finally, integrating out the charge sector yields the ef¬ 
fective action Ss = S^i + S'int in the spin sector at zero 
temperature : 

>S'o=^ 7>s(9,w)(^s(-g,-w) 


1 


^2 + 


_Vs{q,Uj)Ksiq,u}) Ks{q,u}) 


(32) 


5int, = 


2(7ra)^ 


da; dr cos 




The effective spin velocity and the Luttinger parameter 
obey the following equations 


Vs{q,ui)Ks{q,uj) = vp , 


Vsjq.uj) _ Vs f YJL\ 

Ks{q,U!) Ks \2 vfJ 


1 - 






iv*qy 


(33) 
• (34) 


Note that Vs and Ks are the coupling constants of the 
system without SOC. The effective propagator in the spin 
sector is affected in two ways by SOC as can be seen in 
Eq. pi): 


(i) The parameters in the spin sector are explicitly 
renormalized by a Sv term and (ii) there is a contribution 
from the charge degrees of freedom where v* is renormal¬ 
ized according to Eq. (26). 


In the following, we will derive the RG equations of 
the effective action Eq. (32) employing a standard Wilson 


RG procedure in momentum space in order to study how 
SOC affects the phase diagram of interacting electrons. 


B. Renormalization group analysis 


To determine the effect of SOC in the presence of inter¬ 
actions, we study the RG equations of the effective action 
Eq. (32) employing a perturbative Wilson RG procedure 


in momentum space. To this end, we integrate out all 
high energy degrees of freedom between the momentum 
cutoff A and a lower cutoff A' to obtain the low-energy 
physics of the model. Details of the calculation can be 
found in Appendix The result of this procedure is 
encoded in differential equations for the dimensionless 


strength of backscattering yg and forward scattering y 
(the latter is related to Ks). The equations are of the 
BKT type and read as: 


t - 


(35) 


Here, ^ = ln(A/A') and the initial values are given by 


ys{0) = U and y{0) = U + Sv f{U,5v). (36) 


They are determined by the dimensionless strength of in¬ 
teraction U = aoU/TTVp and the dimensionless strength 
of SOC Sv = {Sv/2vf)‘^. The function /(t/, Sv) is nonuni- 
versal in the sense that it depends on the way the inte¬ 
gration procedure of the RG is performed. The universal 
fact, however, is that it vanishes at some point and that 


sign(/(7)) = sign 



(37) 


The ratio of the two dimensionless parameters U and 
Sv therefore completely determines the phase of the in¬ 
teracting electron gas in the presence of SOC. The cor¬ 
responding phase diagram is shown in Fig. In the 
absence of spin-orbit coupling, for Sv = 0, the bare cou¬ 
pling constants are equal, j/s(0) = 2 /( 0 ), and the flow is 
along the separatrix. The SOC changes the initial val¬ 
ues of the coupling strengths and depending on the sign 
of the function f{U,Sv) we either have 2/s(0) > 2/(0) or 
2/s(0) < 2/(0)- For strong SOC, Sv > 2U, the sign of / 
is negative, which results in 2/s(0) > 2/(0) and the sys¬ 
tem will flow to the strong coupling SDW phase. On the 
other hand if Sv < 2U we have 2/s(0) < 2/(0) and the sys¬ 
tem flows to a Luttinger liquid phase with renormalized 
Luttinger parameters. 

To show this, we now review some of the properties 
of the BKT flow equations.!^ The flow equations are in¬ 
variant nnder a sign change of 2/s (but not y) and are 
characterized by the invariant p = yl — y^ and the ra¬ 
tio of initial coupling constants 2/(0)/2/s(0) = cos(/3). We 
point out that while the RG flow is independent of the 
sign of 2/s the nature of the strong coupling fixed point is 
not. Indeed as we will discuss later in Sec. lIII Ci the dom¬ 
inant correlations in the strong coupling phase are either 
of the spin-density wave type (i/s < 0) or the charge- 
density wave type ( 2 /s > 0). 

The flow equations can be integrated to find 

y{l) = HC 0 t{yl +P), 2/s(0 = ■ (38) 

There are three different regimes, illustrated in Fig. 

(I) weak coupling: \y{0)\ > ys{0),y{0) > 0. Here 
l3 = iXi X > 0 and y, = im where m > 0. Backscat¬ 
tering flows to weak coupling, 2 /s 0, while y —>■ m 

as I —/ 00 . This phase is realized for Sv < 2U. 






















SDW (topological) ys 


C. Nature of the strong coupling phase 



CDW (trivial) 


FIG. 3. RG flow of Eq. ( |35[ ) describing the phase diagram 
of interacting ID fermions. In the presence of spin-orbit cou¬ 
pling, the starting point of the flow is changed away from 
the SU(2) invariant line (black rectangle) to the region of 
strong coupling flow (blue circle) for strong spin-orbit cou¬ 
pling, Sv > 2U or to the region of flow towards the Lnttinger 
liquid phase for Sv < 2U (green diamond). 


(II) cross-over Here ^ > 0 and 0 < f3 < tt. The flow 
is still to strong coupling but via an intermediate 
regime (0 < /? < 7r/2) where ys initially decreases. 
This phase is realized for 6v > 2U. 


(Ill) strong coupling: |j/(0)| > ys{0),y{0) < 0. Here 
/3 = ix + 7r,x>0 which yields y, = —im where 
m > 0. Both coupling constants flow to strong 
coupling reaching a pole singularity at (q = x/w. 
This phase is not realized in the context of our weak 
coupling analysis. 


The lines where y = 0 are the SU(2) symmetric lines. 
The flow is to weak coupling for repulsive interaction 
[y(0) > 0] and towards strong coupling for attractive in¬ 
teraction [y(0) < 0]. The RG flow evolves along these 
lines in the absence of SOC. 

We are particularly interested in the study of the 
strong coupling phase where a gap opens in the spin 
sector of the theory. This phase is realized in the case 
of moderately strong SOC, Sv > 2U, i.e., the crossover 
regime in the above classification. 

To get an estimate of the magnitude of the spin-gap 
As, we integrate the RG-flow Eq. (38) up to a scale I* = 


ln(A/As) where gs{l*) = 1. In our case, the deviation 

_2 

from the separatrix is always small and thus y ^ Sv <C 1. 
This yields the estimate 


As = 


(39) 


In this section, we discuss which quasi-long-range cor¬ 
relations govern the behavior of the strong coupling 
phase. According to the Mermin-Wagner theorem, there 
is no spontaneous breaking of a continuous symmetry 
and thus no long-range order in dimensions D < 2 at 
any finite temperature and for short range interactions. 
Therefore, the average of any order parameter {0{x)) 
vanishes. On the other hand, an operator O may exhibit 
quasi-long-range order meaning that the correlation func¬ 
tion ((!I(a:)C>(0)) decays as a power law. The ground state 
of the system is then characterized by the operator whose 
correlations decays the slowest among all. 

We study two possibilities for the order parameter: 

1. Charge-density wave: 

OcDW -I- 

= -[sin(-\/4^(p^) -I- sin(-\/47r(/?^)] (40) 

2 

=-sin(-\/^(/?c) cos{V^fs) ■ 

Tra 

2. z-component of the spin-density wave: 

OSDW + ^IRiJ 

= -^ [sin(-\/47r(/3^) — sin(-\/4^(/?^)] (41) 



Tra 


Superconducting correlations are also possible, but we 
will not consider them here since we discuss the ex¬ 
perimentally relevant case of repulsive electron-electron 
interaction Kc < 1 for which their correlations decay 
faster than density wave correlations. We note, how¬ 
ever, that the crucial element in the following discussion 
is the structure in spin-space (which has two possibili¬ 
ties if gapped), so it is easy to adapt our results for the 
case of attractive interactions, for example, as discussed 
in Ref. [ini 

In the strong coupling phase, the field (ps orders and 
therefore its dual field 9s is totally disordered and all its 
correlation functions decay exponentially to zero (this is 
the case for the x and y components of the SDW and 
we therefore exclude them from our discussion). On the 
other hand, the GDW and SDW order parameters will 
develop quasi-long range order. 

We are now in the position to discuss the nature of 
the strong coupling phase of interacting electrons in the 
presence of SOG. We note that if the cosine term in the 
spin-sector is relevant (massive phase), we have two pos¬ 
sibilities: 


where 0 < /3 < tt. The spin gap is therefore exponentially • < 0, in which case the field Ps wants to lock at 

small. a minimum where cos(-\/^:ps) = I to minimize the 
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energy of the system, i.e. Lps = \/ "k j^n. In this 
case (cos(-\/^(/3s)) 7 ^ 0 and (sin(-\/^(/3s)) = 0, so 
the dominant correlations are of the CDW type. 

• (/s > 0, in which case the field ^ps wants to lock at a 
minimum where cos(-\/^(/Js) = —1 to minimize the 
energy of the system, i.e. ips = \/7r/2(n + 1/2). In 
this case (cos(-\/^(ps)) = 0 and (sin(-\/^(/ 3 s)) 7 ^ 0 , 
so the dominant correlations are of the SDW type. 

In the case of repulsive electron-electron interaction 
discussed above, we have 5s (0) ~ t/ > 0. The system 
therefore develops SDW correlations in the ^-component 
of the spin. 

IV. PROPERTIES OF THE 
SPIN-DENSITY-WAVE PHASE 

A. Effect of disorder 


symmetry of the Hamiltonian constricts the disorder po¬ 
tential to be symmetric in spins, = A. Therefore, 

only the nonmagnetic part of the impurity remains and 
the impurity operator directly couples to the CDW order 
parameter operator: 

^limp = AOcdw • (44) 

Since charge-density-wave correlations vanish in the 

SDW phase, we can already anticipate that impurities 

will be an irrelevant perturbation to the system. This 
statement can be made more precise by the following ar¬ 
gument. In the spin-gap phase, the field pig is locked to 
one of the minima of cos(-\/^:ps)- Fluctuations around 
this ground state can be described semiclassically by ex¬ 
panding 

p)g(x,T)= (p® ■^5p)g(x,T). (45) 


Consider the situation that the system has established 
spin-density-wave order, i.e., the spin sector of the theory 
is gapped. We now want to study the effect of sparsely 
distributed disorder in this phase. This situation can 
be modeled by considering a local scattering potential 
at a: = 0 with strength Ao- for up- and down electrons, 
respectively: 

^ ( Ao-n^ (0) 

G 

= ^ A. [{RlR, + LiL,) + {RlL, + H.c.)] . 

(42) 

The first term describes forward scattering and we will 
neglect it from now on since it is well known that it can 
be removed via a gauge transformation.!^ Thus, 


Ffimp=^A<,(i?^L, + H.c.) 

G 

= -[A^ sin(-\/47r(/3^) -|- A^ sin(V47r(^^)] 

■ sin(-\/^(/5c )cos(v^V^s) (^^) 


7ra 

TTO 

At — A^ • 


TTO 

-^t 4" 


sin(v^V5s) cos{V^(pc) 
OcDW H-• 


In the last line we expressed the impurity operator in 
terms of the CDW and SDW order parameters defined 


in Eqs. (40) and (41). We point out that according to 


the transformation rules in Eq. (13) Ocdw is even under 
time-reversal while Osdw is odd^ote that this does not 
imply that time-reversal symmetry is spontaneously bro¬ 
ken in the SDW phase, since the type of order is not long- 
range). We therefore refer to the two terms as nonmag¬ 
netic and magnetic impurity, respectively. Time-reversal 


In the case of repulsive electron-electron interaction dis¬ 
cussed in this work, the minima are = •y/7r/2(n-|-l/2) 
and thus we get 

■f^imp = sm(V'27r<pc)<5(pg + 0(<5(pg). (46) 

i/7ra 

Integrating out the massive fluctuations S(ps the charge 
sector of the model maps to the Kane-Fisher modeP^ 
with K —>■ 2Kf,. Thus nonmagnetic impurities have no 
effect to first order but can generate a term in second or¬ 
der that can become relevant for very repulsive interac¬ 
tions Kc < 1/2. Magnetic impurities, on the other hand, 
which will only be present when time-reversal symmetry 
is broken, are given by 

.I^imp.magn — COs('\/ 27Vipc) -)- (D(^6(Pg^ , (47) 

TTO 

which becomes relevant for Kc < 2. So even a small 
concentration of magnetic impurities would destroy the 
SDW state. 

The analysis of a single impurity in the SDW phase 
closely resemble s the study of a local perturbation in a 
fermionic ladder.l^SEIl In the latter case, the system con¬ 
sists of two interacting and closely located nanowires with 
CDW correlations which have a relative phase shift of tt. 
The authors find the surprising result that a single impu¬ 
rity placed on a rung of the ladder represents an irrele¬ 
vant perturbation that does not change the conductance 
at zero temperature. On the other hand, two identical 
impurities placed on the same rung represent a relevant 
perturbation driving the ladder to an insulating phase. 
This is completely analogous to the nonmagnetic or mag¬ 
netic impurity in this work. 

Let us now discuss the situation when we are not deal¬ 
ing with single impurities but the impurity density is 
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high. The corresponding Hamiltonian is given by 


a single nonmagnetic impurity the conductance reads as 


Himp — 


E 


dxU{x)na{x) 


G = 


2e^ 

h 


-Ae 


4K^-2 


(51) 


=E 


Uf{x){RlR, + LlL,) 


a 


2 




i 

Itt 


+ {Ub{x)R\La + H.c.) 
dxdx'^cUf{x) 

J dx Ub{x)e~'^‘^'= cos{V¥7r(ps) + H.c.. 


(48) 

Here we dehned the forward and backward scattering po¬ 
tentials 


Uf{x) 

Ub{x) 


f ^ e^<i-u{q ), 

JlgKfcp 


'\q\^kj 


^e^<i-Uiq + 2kF). 

. 


(49) 


where H is a nonuniversal constant and e = max[T, H]. 
To establish the boundary between localized and delocal¬ 
ized regime for the case of weak disorder, it is sufficient 
to replace e —> 1/L and multiply SG by the number of 
impurities Nimp ~ 1/L. Now, if G grows with increasing 
L, the system is in the localized regime or else it is in the 
delocalized regime. This yields the thresshold Kc <3/4 
for localization in the presence of nonmagnetic impuri¬ 
ties and K < 3 for magnetic impurities. This situation is 
very similar to the effect of disorder in quantum-spin-Hall 
edge states.l^ The edge states are also protected against 
nonmagnetic impurities for weak interactions, Kc > 3/8, 
due to the time-reversal symmetry but are expected to lo¬ 
calize if the interactions in the system become too strong. 


B. Luther Emery solution and edge states 


which are white noise correlated Uf{x)Uf{x') = 
Ub{x)U^{x') = D6{x — x'). After averaging over disorder 
the complete model is 


5e = 


2K*v* 


dxdr [{drifc)^ + 




1 


2Ktvt 


Y dxdr [{drip‘^,y + 


9s 

2(7ra)^ 


J dx dr cos , 


SSO =i— Y / <^xdT drifc^x'fl : 

1) T-n ‘ ^ / 


Sf=-^Y J Uf{x)dx(Pc , 


Sb = 


D 

(tto)^ 


Y [ dxdridra cos(v^[(p“(l) - <Pc(2)]) 

a.b 


X cos(-\/27r(/?“(l)) cos(v^</3s(2)). 


(50) 


Here, a, b denote replica indices. First we note that for¬ 
ward scattering can be removed by the gauge transforma¬ 
tion -I- Kcjvc^/n J^dyUf{y). Second, we con¬ 

sider the case where the spin gap already has established 
and we add disorder on top of it. In this case we may ne¬ 
glect the marginal term Sso and expand ips around the 
ground state configuration. Integrating out the massive 
fluctuations Sips the model maps to the Giamarchi-Schulz 
modeP^ with K 2K*. Disorder thus becomes relevant 
for K* < 3/4. 

This result can also be implied from the result for a 
single impurity by the following heuristic reasoning. For 


In the bosonized theory it is possible to investigate the 
properties of the strong coupling phase by refermioniza- 
tion (see Appendix [A|) which enables us to map the sine 
Gordon model at Kg = 1/2 to spinless fermions with 
mass m = ggj2Tra. The Hamiltonian density in the spin 
sector at Kg = 1/2 then reads as 


Rs [n^ + {dxP^f'] + — cos{\/^tp) 

2 Tra 

= - iv{R^dxR - L^dxL) + im{R^L - R) 


(52) 


While this is valid only at one specific value of interac¬ 
tion strength, we will assume that the features of the 
refermionized Hamiltonian characterize the whole mas¬ 
sive sector of the model. One particular property of the 
model Eq. (52) we want to mention at this point is that 


it hosts zero energy modes localized at the edge when 
put on a finite segment 0 < x < L with open boundary 
conditions .132 In the L —>■ oo limit (semi-infinite system), 
the wave function of the zero energy bound state at x = 0 
takes the form 


Xo(x) = -e-^^N (53) 

Vs 

While it is not mentioned explicitly in the original pub¬ 
lication, it is important to realize that the edge state is 
only a normalizable solution if m > 0. 

Since electrons carry both charge and spin, the edge 
state in the spin sector affects the whole electron liquid . 
To see this, consider a boundary between a topologically 
trivial phase with m > 0 (e.g., the vacuum) and the 
topologically nontrivial phase with m < 0 at x = 0. Since 
the field tps is pinned to tpg = \pKj2n\ for m < 0 and to 
(fg = •^7r/2(n2-1-1/2) for m > 0 where ni, n 2 are integers, 
there must be a kink of minimal magnitude a/tt/S in (fg 
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across the boundary. Such a kink in ^ps corresponds to an 
accumulation of half of the electron spin at the boundary: 


Sz = Jdxps{x) = -J= Jdxdx(Ps{x) 


(54) 


where we used the spin-density defined in Eq. (A9). This 
agrees with recent findings in an analogous model in 
Ref. [TUI 


C. Topological classification 


In this section, we show that the the strong coupling 
fixed point in the spin sector of interacting electrons in 
the presence of SOC is a topological insulator of class 
BDI. 

The topological phase of a model of noninteracting 
massive fermions can be determined by the transforma¬ 
tion properties of its single-particle Hamiltonian, which is 
defined in the first Brillouin zone, under the anti-unitary 
operations of time-reversal 0, charge conjugation C, and 
the combined chiral operation 5 ^ 0C. They are defined 
as 

0n*{k)9-^ = + n{-k), e^ = ±i 
Cn*{k)C-^ =-H{-k), C^ = ±l (55) 

En*{k)E-^ =-n{k). 

Consider now the refermionized Hamiltonian ( |52[ ) , which 
describes the effective theory of the spin sector in the 
massive phase. Defining the vector 4' = {R, L)^ the 
Hamiltonian in momentum space takes the form 


H = d'\k){vkaz, — may )m{k) = d'\k)nik)d>{k). (56) 


For the single particle Hamiltonian 'H(fc) in Eq. ( [56| , 
we can explicitly construct the antiunitary operations 
in Eq. (55). They can be represented as 9 = axJC, 


C — K., E = ax, where K, denotes complex conjugation 
and the Pauli matrices act in chiral space. Note, that 
these are symmetries of the refermionized spin sector, 
in particular the operator 9 represents the projection of 
the time-reversal operator onto the gapped spin sector 
and should not be confused with the time-reversal oper¬ 
ator T of the physical electrons defined above Eq. (13). 


From the explicit construction of the operators, we see 
that 9^ = iP = = 1, which places the system into 

the topological class BDI. This class has a Z topological 
invariant in one dimension. 

We now proceed to calculate the topological invariant. 
To this end, we bring the Hamiltonian in the “conven¬ 
tional” chiral form by switching Uz to ax by an appro¬ 
priate unitary rotation. The resulting Hamiltonian reads 
as 


H'ik) 



h{k)\ 
0 ) 


= d{k)a- 


(57) 


where we defined h{k) = vk im and d{k) = 
{vk, —m, 0)"^. We define the normalized matrix element 
of the Hamiltonian q{k) = h{k)/\d(k)\ and the winding 
number 


— / dk q~^ {k)dkq{k) 
J — OO 


-sign(m) . 


(58) 


The winding number is not an integer since it is origi¬ 
nally defined as a mapping from the compactified Bril¬ 
louin zone, isomorphic to 5^, to Z. The low-energy the¬ 
ory defined on a non compactified manifold (R) does 
not capture the physics in the whole Brillouin zone but 
only at one time-reversal-invariant point [the other half 
of the winding number is contributed by the other time- 
reversal-invariant point but we have no means to tell if 
it contributes sign(m) or —sign(m)]. From the winding 
number we can therefore not deduce which sign of the 
mass corresponds to the topological and which to the 
trivial phase. We can, however, determine the number of 
edge states at the boundary of two topological materials 
with masses mi and m 2 . It is given by n = v{mi) — i'{m 2 ) 
and since i/ 0, we know that one sector has to be topo¬ 
logically nontrivial. In fact, we have already established 
that the SDW phase with m < 0 hosts zero energy edge 
states and therefore this has to be the topologically non¬ 
trivial one while the CDW phase is a topologically trivial 
insulating phase. Note that the Z topological invariant 
can take only two values in the present case because we 
consider only a single quantum wire which can either 
have edge states or not. If we considered multiple identi¬ 
cal wires, scattering between them cannot gap the indi¬ 
vidual edge states as long as the anti-unitary symmetries 
are preserved leading to the Z topological classification. 
This is similar to the topological properties of multiple 
copies of the Su-Schrieffer-Heeger modeP^ for polyacety¬ 
lene. We want to stress that the above classification 
scheme applies only to the gapped spin sector, which can 
be mapped to noninteracting fermions. Nonetheless, it 
allows us to clearly identify that a topological and a triv¬ 
ial phase exist. The properties exhibited by the whole 
electron liquid if the spin sector is in the topological phase 
have been discussed in the previous sections. 


V. SUMMARY 

We investigated the influence of spin-orbit coupling 
in interacting one-dimensional quantum wires. Spin- 
orbital interactions break the spin-rotational symmetry 
from SU(2) to U(l) and break the inversion symmetry in 
the wire. This lifts the spin degeneracy of the spin-up and 
-down energy bands and leads to different Fermi veloci¬ 
ties vi ^ V 2 in each band which is depicted in Figj^ In a 
specific hopping model discussed in Sec. [Hj we estimate 
that the velocity difference is of the order of 5v ~ at' jt'^, 
where a is the strength of SOC and t, t' are the am¬ 
plitudes of nearest- and next-nearest neighbor hoppings. 
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respectively. We find that this difference in velocities 
drastically affects the nature of interacting electrons. 

The most interesting situation occurs when the dimen¬ 
sionless strength of SOC is stronger than the strength of 
interactions 6v > 2U. In this case, interactions drive 
the spin sector of the system to a strong coupling phase 
with quasi-long-range spin-density-wave order where a 
spectral gap is dynamically generated. This prediction 
was established by treating interactions using bosoniza- 
tion and a weak coupling renormalization group analysis 
which is controlled by the small parameters of dimension¬ 
less interaction U and spin-orbit-coupling strength 6v. 

We show that the spin sector of the gapped spin- 
density-wave state is topologically nontrivial (symmetry 
class BDI). Meanwhile, the charge sector of the quantum 
wire is still massless. Since physical electrons carry both 
spin and charge the whole electron liquid forms an un¬ 
usual topological state. The topological nature manifests 
itself in the emergence of zero-energy edge modes at the 
boundary of a finite system that carries fractional elec¬ 
tron spin. The existence of these modes is protected by 
the bulk spin gap. Furthermore, we find that due to the 
gap in the spin sector, the bulk of the system is protected 
against nonmagnetic impurity scattering for interaction 
strengths Kc >3/4, for high impurity concentration and 
for Kc > 1/2 in the case of a single impurity. These 
extraordinary transport properties are protected by the 
time-reversal symmetry in the system. 

We hope that our work will stimulate experimen¬ 
tal search and studies of strongly interacting one¬ 
dimensional systems with sizable spin-orbit coupling. In 
particular, the robustness against disorder and the gap¬ 
less edge states of the system should be probable in ex¬ 
periment. We also note that this robustness to disorder 
may eventually be technologically exploitable in nano¬ 
electronics, as the near-perfect conduction property of 
the wire is protected by topology and does not require 
ultra-clean wires. 
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Appendix A: Bosonization conventions 


We use the following conventions: 


^ p-»V47r<A^(a:) 


Raix) = . - e 

vStto 


(Al) 


where a is the bosonic UV cutoff. The Klein factors an¬ 
ticommute {ko-jKo-’} = 26cr,a-' and ensure correct anti¬ 
commutation relations of fermions with different spins. 
Since they are not dynamic variables, we are free to 
choose a specific representation: 

kI = 1, K^Ki = -KiK^ = i. (A2) 

The correct anti-commutation relations of right- and left- 
movers are fixed by the commutations relations between 
(j)^ and (j)^: 




(A3) 


=-'n6^^n'^a,o' sign{x - y). (A4) 


and the Campbell-Baker-Hausdorff rule 




(A5) 


which holds when the commutator of A and B is not an 
operator. We also introduce the conjugate variables 

a — A.L aR 




(A6) 


where 11 = dxd is the conjugate momentum to tp. Using 
the above definitions, we find 




r],o 

- 


Ul T — 

1 


(A7) 


RlRcr + LILct = -^dxPa ■ 

VTT 


The spin and charge degrees of freedom are given by 


+ Vi Vt- Vi 


(AS) 


In terms of these operators, the nonoscillatory part of the 
charge density and the z-component of the spin density 
read as 


Pcix) = ^'ilji{x)il}c{x) = —^dxVc{x), 

(7 ^ 

pI{x) = '^'4’l{x)Slcr'i^^,{x) = ]^'^aii)l{x)^^{x) 

(7 <7 

= , _ dxVsix) ■ 

(A9) 


Finally, we state the following refermionization identities 
for spinless fermions: 



0 —i2naL^R, 

(AlO) 

giV 47t9 

O i2Tra R^ , 

(All) 

dxV 


(A12) 

n 

0 - A^'A) , 

(A13) 
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Appendix B: Details of the momentum space RG 


We introduce the vector notation q = {q,ui/vs) and 
r = {x,Vst). The absolute value is denoted by r = |r|. 
To perform the RG analysis of Eq. (32) we split the fields 
into fast (>) and slow (<) components (we drop the sub¬ 
script s from now on since we exclusively deal with the 
spin sector): 


ip{x,T) = ip>{x,T) + ip<{x,T) , (Bl) 

where we defined 

= I ^ (B2) 

|q|<A/fc 

ip>{x,T)= j ip{q,w) . (B3) 

A/&<|q|<A 


Now, we integrate out the fast degrees of freedom whose 
momenta lie between the momentum cutoff A and the 
new smaller cutoff A/&, defined by the scaling factor b = 
« 1 -I- £ > 1. This yields the effective action 


=S^ + (5int)> - ^ 




(B4) 


Sint = 


2(7ra)' 


Jdxdr cos 


(B5) 


To perform the average, we need the correlation function 
of spin fields in momentum and frequency space which 
can be obtained from Eq. (32) and reads as 


K, 


V. 


's {u)/vsY + q^ - q^5 


1 - 




(u;/Rs)2+'y2g2 

Here, we defined the dimensionless constants 


(B 6 ) 


7 = — , and 

Vs 




(B7) 


The velocities in the charge and spin sectors are defined 
in the main text in Eqs. (20) and (26). They can be 
expressed in terms of microscopic parameters as 



ttvf ’ 


and 



(B 8 ) 


In order to make analytical progress, we assume weak 
interactions U 1 and weak SOC du <C 1. This allows 
us to make the expansion 


7 - l + U- -dv, 


and (5 ~ . 


Eurthermore, we define the functions: 


(B9) 


Fj.si^r) = — {(p>{x,T)(p>{0,0))^ , 
A^^siAr) 

^ g-2<?Ar,F-,.i(Ar) 


(BIO) 

(Bll) 


In general, both F^ siv) and A^^s(r) will depend on 
the parameters 7 and 6. Before we proceed to derive 
the RG equations we shortly comment on the properties 
of the function F^^siv). We first consider the case of 
(5 = 0. Then, we find F{Ar) = Jo{Ar). The zeroth 
Bessel function of the first kind is oscillating and falls 
off as a power law as r —)■ cx). This is, however, not the 
behavior we would like to have since it does not allow us 
to perform a gradient expansion in r. It was shown that 
this is a consequence of the fact that we chose a hard 
cutoff for the radial integration. If we were instead to 
choose a smooth cutoff F{r) is truly short range.l^In the 
following, we will not specify the explicit form of F{r) 
but just assume that it falls off rapidly as r —> 00 . 


We now calculate the funct ion F a s{0) which turns out 
to determine the flow in Eq. (B31): 




dqduj Kg 


Vs q2 — q2§ 

1 


A/6<|q|<A 

I ^ _ 

^ J A/b 9 Jo 27r 1 —(5cos^(0) 

d0 1 

27T 


1 - 


4{io/v,y 


(a;/Rs) 2 + 72 g 2 


1 - 


4sin^(6)) 


sin^(0)+72 cos2(0) 


/o I —i5cos^(0) 


1 - 


1 + 7 ^ cot‘^{ 6 ) 


(B12) 


Analytical solution in the limit of small 6 


The integral in Eq. (B12) cannot be performed analyt- 
































14 


ically but we can find a analytical expression in the case 
of 5 ^ 1. To first order in 6 we find: 


r'^^de 

./n 27 r 




1 - 


1 + (5 cos^( 0 ) 
4 


1 - 


(1 + 7 )^ 

In the last equation we defined 


1 + 7 ^ cot^( 0 ) 
= l + ^f(7)- 


(B13) 


f(j) = ^-^-. 

2 (1 + 7)2 


(B14) 


Notice that f(-j) > 0 for 7 > 1, f(j) < 0 for 7 < 1 
and /(I) =0 (at this value of 7 the lowest non vanishing 
correction is of order 6 ^). This can be summarized as 


sign[/( 7 )] =sign 



(B15) 


Derivation of the RG equations 


We now proceed to derive the RG e quations by calcu¬ 
lating the effective action in Eq. (B5). To first order we 
find 


^2(f^ yds dr (cos [y8)((v3<+ +>)])^ 

= —r~T^ /da;dr cos(-\/87r(/?^) 

2(7ra)^ J 

= (2 - 2KsF^^siO))i [dxdr cos{V^ip). 

2 {Tray J 

(B16) 


In the last step we rescaled space time (a;, r) —)■ b{x,T), 
so that the theory is again defined with the cutoff A. We 
now define the short hand notation 1 = (a:i,ri). The 
second order is 



^4(la)4 (+^++^)(l)] cos [v^((p<+(p>) (2)])^ 

- (cos [v^ (+< + (1)] )> (cos [V^ {if< + ip>) (2)] I 

l_3s_ /dld2 ^ giN/^(«iV^(l)+e 2 V^( 2 ))| ^gi%/+f[£i¥;>(l)-|-e 2 <,= ^( 2 )]^^ _ g-87r([¥>>(0)]^)> 

2 (27ra) 

(2^^^)4 ^~^(0) ydld2 |cos (v^[(/j<(1) - (/j<( 2)]) [^-^(Ars) - l] 

+ cos (^v^[(^<(l) + +^( 2 )]^ [A 2 (Ars) - l] } . 


(B17) 


where the function A(Ar) is defined in Eq. (Bll). We 
now introduce relative and center coordinates 


AT '2 2; = a;i-a;2, (B18) 

r = ^ , t = ti-T2 (B19) 


We expand the terms in brackets to first order in i 

A-2(0) [A-2(Ar) - 1] ~ 2«,E„,5(Ar«), (B21) 

A-2(0) [A2(Ar«) - 1] ~ - 2£K,F^^s{Kr). (B22) 

Since F{r) is decaying rapidly as r —>■ 00 we can perform 
a gradient expansion of the cosine terms. 


and the vectors R = {X,VsT) and r = {x,Vst). This 
gives 


cos ^•\/^[:/j^(R + r/2) — (/ 5 ^(R — r/2)]^ 
~ 1 - 47r(^VR(/5(R)r^ , 


(B23) 


S. 


( 2 ) 

eff 


I cos ( 
+ cos 


{2TraY 


rA-2(0) d'^Rd^r 


'^[(/j‘'(R + r/2) — (/?^(R — r/2)]^ [A ^(Ar) — l] 

/^[(^“^(R + r/2)) + (/j<(R - r/2)]^ [^^(Ar) - l] | 

(B20) 


cos 


/^[(/3^(R + r/2) + +^(R 
cos ^•\/327r(/5^(R)j . 



(B24) 


The cosine term that appears is less relevant than the 
"original cosine term and we therefore neglect it. The 
first term however gives the following correction to the 
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quadratic action 


( 2 ) ^ 

(27ra)4i;2 


£ Jd^R(vRipi-R)y Jd\ r^F^^siAr). 

(B25) 


The integral over relative coordinates gives a nonuniver- 
sal number Ai. We define 


It is straightforward to show that, to first order in g^, the 
RG equations for yg are not affected under this redefini¬ 
tion. The total set of flow equations then reads as 


= [2-2KMF^,siO)]ysi£), 

dK^je) ^ ys{£?Ks{£)^ 

d£ 2 


(B31) 


/ 27r 

d\ r^F,4Ar) = -^A,{^,6). (B26) 


Since the correction is already linear in £ we can set the 
rescaling factor of the coordinates b to one. The resulting 
correction reads as 


c(2) _ 


g'^sKs 


{aA)^'K‘^v, 


■Ai{'^,S)£ dxdr [dx^p) 


{drpf . 

(B27) 


The term gives a correction to Vs/Kg but leaves l/vgKs 
invariant, which means that only Kg will flow. From 
Eqs. (B16) and (B27[) we find the differential equations 


dKg 

d£ 

dgg 

d£ 


gl{f)KK(-) 

[aAYvl 


^ 1 ( 7 ,^) 


{2-2Kg{£)F^4£)))gg{£). 


(B28) 

(B29) 


We can bring them in the conventional BKT form by 
defining the dimensionless coupling constant 


y 


2 

S 


^glKg 

{aK)'^TT‘^vl 


Ai(7,(5) . 


(B30) 


Notice that these equations are perturbative in gg but 
exact in Kg and the parameters 7 = v*^/vs and 5 = 
{dv/2vg)^. 

We simplify the flow equations for small interaction 
strength and small S by setting Kg{£) = 1-1- g{£)j 2 and 
using the asymptotic form of F-^^s derived in Eq. (B13|. 
This yields 


dyg 

d£ 

dy 

d£ 


- ys{£){g{£) + 

-yl{^)- 


d /( t )) + 0{ygg5) , 


(B32) 


We can define y{£) = g{£) -I- 6 /(y) which leaves us with 
the conventional KT flow equations. 

^ =-yg{£)y{£), (B33) 

f = -yK^)- (B34) 

The initial values of these equations are given by yg{0) = 
U and y{0) = U + Sv Note that in the absence 

of SOC {S 0) the RG equations have to coincide to 
describe the flow along the separatrix of the BKT flow. 
Therefore, we must demand that g{0) = yg{0). 
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